Universality of the one dimensional Bose gas with delta interaction 
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We consider several models of interacting bosons in a one dimensional lattice. Some of them are 
not integrable like the Bose-Hubbard others are integrable. At low density all of these models can 
be described by the Bose gas with delta interaction. The lattice corrections corresponding to the 
different models are contrasted. 



o 

CD 



03 



i 

C 

O 
o 



CO 

> 
in 
in 

in 
o 
^i- 
o 

-I— » 



i 

-a 
c 

o 
o 



X 



I. INTRODUCTION 

The Bose gas with delta interaction is one of the 
most famous models for strongly interacting bosons Q, @ 
(the corresponding classical dynamics is described by 
the non-linear Schrodinger equation). In one dimension 
this model possesses infinitely many, mutually commut- 
ing integrals of the motion and is therefore called inte- 
grable. The integrability of the Bose gas can be for- 
mulated within the Quantum Inverse Scattering Method 
(QISM) that exploit a crucial property of the two- 
body scattering matrix of the system: it satisfy the Yang- 
Baxter equation. This reduces the many-body dynamics 
to a two body dynamics. In the Appendix we explain 
the QISM & the algebraic Bethe Ansatz [it uses a trans- 
fer matrix]. In the last decades many authors studied 
integrable lattice regularization of the continuous Bose 
gas llllll, 00 HE [H El El El The procedure it 
is not straightforward since a direct quantization of lattice 
classical integrable theory whose continuous limit is the 
classical continuous Bose gas field theory leads to "non- 
local" Hamiltonians (that is the interaction involves far 
neighbors). This inconvenience was circumvented in two, 
somewhat complementary ways: keep the scattering ma- 
trix of the continuous theory unaltered and change the 
transfer matrix and then the Hamiltonian 0, HH 0j 
alternatively, change the scattering matrix keeping the 
Hamiltonian formally unaltered H E|. In both cases, 
however, the Hamiltonians though (quasi) local, involve 
higher order processes whose physics is hard to read-out. 

Another paradigmatic lattice model for interacting 
bosons is the Bose-Hubbard model. It was proposed 
first by Haldane E3 ( see a l so E3)> an d widely used 
in many different contexts. In mesososcopic physics it 
describes granular superconductors, Helium films, and 
Josephson junctions arrays |l7j . In quantum optics, the 
Bose-Hubbard dynamics was proposed to model bosonic 
atoms in optical lattices E3 and as effective theory for 
higher spin chains E3l- The Bose-Hubbard Hamiltonian 
describes charged bosons hopping in a ci-dimensional lat- 
tice and experiencing a Coulomb repulsion. The resulting 
competition between kinetic and electrostatic energies in- 



duces quantum fluctuations that dominate the zero tem- 
perature phase diagram. At commensurate number of 
bosons per site the Bose-Hubbard model predicts a Mott 
insulator-superfluid quantum phase transition controlled 
by the ratio between the bosonic tunneling rate t and 
the Coulomb repulsion. In optical applications the sys- 
tem parameters are tuned by changing the laser intensity 
|2fj|: the criticality of the superfluid- insulator quantum 
phase transition is substantially altered by the confining 
harmonic potential |2l| . 

In this paper we focus on one dimension. We remark 
that d = 1 is more than an academic case for the physi- 
cal systems captured by the models under consideration. 
In fact bosonic chains are commonly realized both with 
optical confinement [2(1 ] and, using a different technol- 
ogy, with Josephson junctions E3- Quantum effects are 
strongest in the one-dimensional case. For example, in 
the dilute limit of optical lattices quantum many body 
effects cause the breakdown of the Gross-Pitaevskii the- 
ory (valid for 3d systems) and drive the system to the 
so called "Tonks- Girardeau" regime |22| where the exact 
results of Lieb and Liniger are crucial 23] . Exact results 
of the Bose-Hubbard model would be desirable to study 
physical systems like the above ones, at higher densi- 
ties. However the exact solution of the one dimensional 
Bose-Hubbard model is still a challenging problem. For 
infinite range hopping the thermodynamics is fully estab- 
lished [25j . For nearest neighbor hopping, the coordinate 
Bethe ansatz fails [24| due to the multi-occupancy char- 
acteristic of the bosonic statistics |26j . Recent numerical 
studies suggested that the distribution of level spacings 
obeys Wigner-Dyson statistics |27l , this considered as in- 
dication of non- integrability 28] . 

In this paper we discuss the integrable continuous 
field theory of the Bose gas as universal low density 
limit of various integrable (Faddeev-Takhtadjan-Tarasov 
B. Izer gin- Korepin and quantum Ablowitz-Ladik 
SIlEllllEi) and non integrable (Bose-Hubbard) lat- 
tice models. In particular a number of results obtained 
for the Bose gas can be used for the behavior of the sys- 
tems phrased by Bose-Hubbard models. Finally we dis- 
cuss the differences between the existing integrable lat- 
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tice regularization of the Bose gas and the non-integrable 
Bose-Hubbard model. 

The paper is laid out as follow. In the next section we 
summarize a number of properties of the integrable field 
theory of the Bose gas. In section III we apply some of 
these properties to the Bose-Hubbard model. The lattice 
integrable theories are discussed in section IV. In section 
V we draw our conclusions. In the Appendix we summa- 
rize the main ideas of the QISM. 



II. THE BOSE GAS WITH DELTA 
INTERACTION 

We shall start with continuous field theory. This will 
be the continuous limit of all bosonic models considered 
in this paper. The one dimensional Bose gas is described 
by an integrable non-relativistic model of field theory 
3]. The Hamiltonian is: 

H BG = J dx[{d x ^){d x y) + ctfttf 1 "** - htfttf] . 

(1) 

Here: [*(x), #%)] = S(x - y) and [&(x),&(y)] = 0. 
The coupling constant is denoted by c, and h is the chem- 
ical potential. It is equivalent to a many body quantum 
mechanical problem with Hamiltonian: 



£)2 

H ^ = -J2d^ +2c E 8{ Zj -z k ). (2) 

j=l 3 N>j>k>l 

The energy is E = Ylj=i ( u | — h) > tne momenta u's obey 
Bcthe equations: 



jL _ yj uj -u k + ic 



n 



Uj — Uk — ic 



j = l...N 



(3) 



Let us consider repulsive case c > 0. In the thermody- 
namic limit N, L — > oo such that N / L = D = const the 
ground state energy-density [per unit length] is 

/<? rq 
e (u)p(u)du = / (u 2 — h)p(u)du . 
-q J —q 

The function p(u) describes the distribution of momenta 
in the ground state. It is defined by the equation: 

2irp(u) = 1+1 K(u — v)p(y)dv == Z(u) , 

J-q 

where Z(u) is called dressed charge. The kernel of the 
integral equation is K(x) = 2c/(c 2 +x 2 ). For further pur- 
poses it is convenient to define a function t(u) through 

r 

2ttt(u) = K(u — q) — dvK(u — v)t(v) . 

J-q 



The density, given by D = J" p{u)du, depends on the 
chemical potential since q = q(h) 0- The Bose gas is 
characterized by different phases depending on the sign 
of the coupling constant c. Let us describe the ground 
state at zero temperature: 

i) For c = we have free bosons. In the ground state all 
the particles have zero momentum (Bose-condensation) . 

ii) In the repulsive case c > the Pauli principle is 
valid (this is the distinctive feature of the one dimensional 
bosons). The ground state is a Fermi sphere. 

Hi) In the attractive case c < the ground state is a 
single large bound state of all the particles. All the par- 
ticles stay close to one another (in configuration space). 
It is a droplet. One can say that in this case the phases 
separate. 

In the repulsive case at zero temperature the asymp- 
totics (large distances) of correlation functions is |3j, |29| 

(*t(a;)tf(0)) = Ax- 1/e (4) 

(*t(z)*(:z)*t(o)*(0)) = D 2 + 4 + B COS[2 l Dx] , (5) 

x z x" 

where the critical exponent 9 depends on the density in 
a remarkably simple way through: 

9 = 2Z 2 (q) = 4ttD/v f . 

The thermal correlation functions can be obtained by re- 
placing: x — > ^ sinh (ttTx/vf) in these expression. Here 
vp is the Fermi velocity [for the model under considera- 
tion, it coincides with the velocity of sound]. So in the 
repulsive case the asymptotic of the thermal correlation 
functions is 

(*t(z)*(0)) T - (H) ^ exp (-ttTx/v f 9) , (6) 



(* t (x)*(x)* t (0)^(0)) T = B 2 



2nT 

v F 



^-2ttTx/vf 



2ttT 

Up 



D -irT9x/v F 



cos(2ttDx) . (7) 



The coefficients A, B, B2, B% are related to certain form 
factors 1291. 



III. APPLICATION TO THE BOSE-HUBBARD 
MODELS 

The Hamiltonian for the one dimensional Bose- 
Hubbard model reads 

H B h = ^2 [U(rii-1)- p]ni-t(alai +1 + al +1 ai) , (8) 

i=-N s 

where the operators rii := a] a, count the number of 
bosons at the site i; operators a,-, aj obey the canoni- 
cal commutation relations [of, aj] = <5y and and 2N S + 1 
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is the number of sites. The parameters t, U of iJSJ are 
the hopping amplitude and the strength of the on-site 
Coulomb repulsion, respectively, while the chemical po- 
tential fj, fixes the average number of bosons in each 
site. The density of bosons in the lattice is given by 
D = N/(2N S + 1)A. Here A is the lattice spacing. The 
Bose-Hubbard model can be considered as a possible lat- 
tice regularization of the Bose gas. In fact at small filling 
factor v = N/(2N S + 1) = DA the Bose-Hubbard model 
can be described by the Bose gas with delta interaction. 
In order to obtain the correct commutators, in the con- 
tinuous limit A — > one has to re-normalize the opera- 
tors: <n = y/A®(x) and n t = A^(x)^(x), x = Ai. The 
Bose-Hubbard model then reduces to the Bose gas 



are interesting by them-self: good examples are Calogero- 
Sutherland, Haldane-Shastry [3(J, and pairing [3lll3a . l3^ | 
models. To obtain local Hamiltonians the Lax matrices 
have to be modified. This modification induces complica- 
tions on the final form of the though local Hamiltonians 
(see the Refs. 0, IS 0.0)- All these are solvable by alge- 
braic Bethe Ansatz (see Ref. Q); some basic formulas 
of the QISM are summarized in Appendix. Here we will 
present these Hamiltonians. Some of them simplify in the 
weak coupling limit. In this limit the resulting Hamilto- 
nians are compared with the Bose-Hubbard model. 



Izergin-Korepin models 



Hbh — tA Hbg j 

with c = U/{tA) and h = -Qj, + 2t)/(tA 2 ). Therefore 
results for the Bose gas can be phrased for systems cap- 
tured by the Bose-Hubbard model. The asymptotics of 
the zero-temperature correlations in the repulsive case 
are 



1/0 



(9) 



, v A cos \2TTv(i — j)/A] , . 

a {1-3) {i-jr 



where the critical exponent 9 is 



2 1 



AD 



2 1 



4vt 
IT 



The thermal correlation functions read 
( a l a j)r = (^) ' sinh[7rT(i- j)/v F ]~ 1/6 



+ B 3 



V V F J 

2nT" 



(11) 

(12) 
(13) 



^T0(i-j)/ Vp [ 2 TTv(i-j)/A] 



IV. INTEGRABLE LATTICE MODELS 

In this section we discuss the integrable bosonic theo- 
ries arising as quantum lattice regularization of the field 
theory of the Bose gas with delta interaction. From 
a classical side the Hamiltonian structure [action-angle 
variables] is directly related to the i?-matrix. This con- 
stituted motivation to define the lattice models keeping 
unaltered the i?-matrix (for an alternative way to con- 
struct lattice regularization of the Bose gas see the Ap- 
pendix 0,0)- Pursuing a quantum version of this proce- 
dure, namely directly quantizing the classical Lax matri- 
ces, it turns out that the arising Hamiltonians are non- 
local. Here we should note that non-local Hamiltonians 



The first integrable version of the Bose gas on the lat- 
tice was constructed in Refs jj, |S| ■ To write the Hamil- 
tonian of the model in a compact way we introduce the 
following notations. The bosons interact differently in 
odd and even lattice sites. Let us start with canonical 
Bose operators 

[a j ,a{] = 6i ■ a i |Q) = 0. (14) 

In order to describe the model, it is convenient to intro- 
duce 'renormalized' operators 



bj = a 3 <j j 



(15) 

Here tjj is slightly different in odd and even lattice sites: 



if j is even , 



| yd + aj(ajaj — 1) if j is odd . 



(16) 



where u> = cA/4. We can express everything is terms of 
'renormalized' operators bf. 



(17) 



a) = {\-ub\bj}- 1 for j even, (18) 

For odd j the expressions are slightly different: 

a]a 3 = (1 - w)b]bj{l - wtf bj}- 1 , (19) 

a) = (1 - w){l - wbtbj}- 1 for j odd . (20) 

We can also write a closed commutation relations for 
renormalized operators: 

1 if j is even 

(1 — u))^ 1 if j is odd . 



(21) 



In different lattice sites operators bj commute. Now we 
can present the Hamiltonian of the model as 



H 



IK 



Aa-aj 




(22) 
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Below, we shall represent operators tj in the form 

1 



U 



(23) 



The formula for operators O and C are different for odd 
and even lattice sites. For j odd: 

Oj = <Tj+x{\- + wb\ +l bj+2}(Tj+2, (24) 

^Vj-A 1 - w &]-i & i+i}oj+i > ■ ( 25 ) 



For j even: 



Oj 

c, 



cTj_ 2 {1 + w6)_ 2 6j-i}o-j_i , (26) 



{1 - ub]_ 1 b j+1 }a j+1 . 



(27) 



In the weak coupling limit the Hamiltonian describes a 
coupling of five neighbors lattice sites j — 2 . . . j + 2. Re- 
taining contributions up to second order in cA it has the 
form 



Hik — 



3cA 3 



(28) 



hj = ^2 \vat3n J+a rij + p + w a pa 3+a a^ +l3 (29) 



where 2£j m = ctja^ + a m a, (with certain coefficients 
Va/3 j , i aj g j 9q/3)- After simple manipulations 
Eq. I|28fl can be written as 



cA cA 
= ^Bff - 9—K^ w + g( — ) 2 [K^_j_ + hj] , 



(30) 



where periodic boundary conditions have been consid- 
ered: 



J2 nj = Yl n i+i and 53 n ? = 53 



(31) 



the parameters in Hbh above are fi/g = (cA)(l — 
5cA/16)/8, U/g = 5(cA) 2 /128, t/fl = cA / 8 with 
g = — 4/(3cA 3 ). The Bose-Hubbard model differs from 
Itzergin-Korepin model in the presence of the quadratic 
hopping and in the non local terms. In the continuous 
limit, when the lattice spacing A — > 0, this lattice model 
turns into a continuous Bose gas JTJ. 



Faddeev-Takhtadjan-Tarasov models 

In this section we discuss another lattice integrable ver- 
sion of the Bose gas. The Hamiltonian was suggested by 
Faddeev-Takhtadjan-Tarasov in Ref. Q. The model has 
the same i?-matrix as the continuous Bose gas. It looks 
like a quantum spin chain with negative spin. Let us first 
introduce an operator of 'angular-momentum' Jj,j+i- It 
is defined as a solution of the following operator equation: 

J jd+1 (Jjj + l + 1) = 2Sj <g) 5j+i + 2s(s + 1). 

Here the 'spin' is s — — 2/(cA). The Hamiltonian of the 
Faddeev-Takhtadjan-Tarasov model is 

h ftt = -2 K ^r'(j 7 - i+1 + i)/r(j JJ+1 + 1) . 



The relation with the lattice bosons becomes transparent 
using the Holstein-Primakov realization of spins: 



S? = 



*7 



l ',Pi ■ f>,"j 

V^A ' 
-atpj+pjCLj 



(32) 



57 



cAa]a 



3^3 



cA 



with p = y/l + cAata/4. For small c, the spin s is large. 
The Hamiltonian simplifies (second order in 1/s): 



H 



FTT 



-^[(aj) 2 + o] +1 a] + (a] +1 ) 2 ][a 2 + a 3 +ia 3 + , 
+^a] +1 a]a j+1 aj - s 2 a]a j . 



'j+U 



(33) 



The model is solvable by algebraic Bethe Ansatz |7J . Af- 
ter simple manipulations Ea. (|33J) can be written as 



Hftt = Hbh 



(34) 



n 3 n 3 + l + ~ K j,j+l + K + n i+l) K j,j+l 



where periodic boundary conditions have been consid- 
ered (see Eqs.Jjnj); the parameters in Hbh above are 
fi = k(cA)(cA/16 + l)/4, U = k(cA) 2 /64, t = kcA/8, 
and V = k(cA) 2 /64. Eq. elucidates the rela- 

tion between the Faddeev-Takhtadjan-Tarasov and Bose- 
Hubbard models. In the continuous limit, when lattice 
spacing A — * the Faddeev-Takhtadjan-Tarasov (and 
the Bose-Hubbard model as well) lattice model turns into 
the continuous Bose gas (Q. 
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V. CONCLUSIONS. 

We considered several models of interacting bosons. 
Faddeev-Takhtadjan-Tarasov model and Izergin-Korepin 
model are different lattice models originally constructed 
from the Bosc gas continuous field theory. For these mod- 
els the lattice contributions do not destroy the integrabil- 
ity of the original field theory; instead the Bose Hubbard 
model contains lattice-effects leading to a non integrable 
dynamics. In the weak coupling limit these integrable 
Hamiltonians and the Bose-Hubbard Hamiltonian can be 
compared explicitly. In particular the model in Ea. H34|) 
could serve exact studies of chains of Josephson junctions 
away from the degeneracy point 17]. The lattice effects 
vanish at low density: all these theories are integrable 
and described by the Bose gas with delta interaction. In 
this regime the known asymptotic of the Bose gas correla- 
tion functions can be used for physical systems captured 
by the Bose-Hubbard model (see [38( for a recent applica- 
tion). The Luttinger liquid field theory in Ref. [3j| and 
the Bose gas are characterized by the same set of correla- 
tion functions where the spectral parameters vj,vf,vn 
are determined by (3?} 9 = 2y / vj/vpf, vf = y/vtfVj, and 

vp = e' (q) — f_ due' (u)T(u) /[Z{q)\ (see section II). 



Discussions with R. 
are acknowledged. 
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is the analog of the Casimir operator of Lie algebras and 
it is another invariant of the Yang-Baxter algebra. In 
Eq. IpTTjl V is the projector 



V = R( mV ) :-- 



1 



2 cosh 7/ 



/ 

e~i 

-1 

\ 



-1 
e 1 ? 
0/ 



(38) 



To construct lattice integrable theories for the Bosc 
gas two different procedures (somehow complementary) 
have been pursued. The first consists in keeping the R- 
matrix of the continuous theory unaltered and changing 
the Lax operators. Such procedure has been followed in 
the text above. One can also change both the R and 
L operators to simplify the final form of the Hamilto- 
nian. This approach was initiated by P.P. Kulish et al. 
0,0 and developed in the papers 0,E3>0]> quantizing 
classical discretization on non-linear Schroedin ger equa- 
tion developed earlier by Ablowitz and Ladik [35|. For 
the Ablowitz-Ladik model the Lax operators and the R- 
matrix MM are : 



HO 



Here exp c 



-4 r 1 



/ 1 o o o \ 

b~ c 

c b+ 

\ 1 / 



(39) 



£/£ and £,£ are spectral parameters, 
e ±7? sinh d> 



APPENDIX 



sinh [(p — rj) 



The construction of integrable lattice models of inter- 
acting bosons is based on Quantum Inverse Scattering 
Method (QISM) 0. Here we shall summarize some fea- 
tures of this method. The starting point of the QISM is 
a local quantum Lax operators Li(u) and a matrix R(u), 
satisfying the Yang Baxter equation 

R(<t>)HQ <8> HO = Ht) ® HC)R(4>) ■ 



The monodromy matrix 

T(C) = L_ Ns (C)i_jv a +i(C) • ■ ■ L Ns (C) , 
also fulfills the Yang-Baxter relation 



(35) 



R(<f>)T{C) ® T(0 = T(0 ® T(QR(<P) . (36) 

The transfer matrix is defined as := tr^T(^) where 
tr(o) means the trace in the auxiliary space. It is a gener- 
ating functional of integrals of motion and of the Hamil- 
tonian since it commutes with itself at different values 
of spectral parameters: [i(C)>*(£)] = (*(C) is an invari- 
ant of the Yang-Baxter algebra). The quantum determi- 
nant |2J| 



det, (T(C) :=TrVT(C)®T(e-X) 



(37) 



sinh rj 



sinh ( 



v) ' 



the quantity rj € IR is a deformation parameter. The 
Hamiltonian of quantum Ablowitz-Ladik can be defined 
by 



- E 

j=-N s 



qj+i + q) +1 qj - aiog(i + qW, 



(40) 

We shall see later that a doesn't play the role of coupling 
constant. In fact the potential and the kinetic energies 
commute (then the potential energy turns into a con- 
stant on the eigenstates of the kinetic term of (|40|) ). The 
operators qj satisfy the following commutation relation: 



[q l ,q}] = (e 2 ^l)(l + q}q J )d 



(41) 



as follows the Yang-Baxter algebra. The Hamilto- 
nian (|40() can be rewritten using the trace identities. In 
particular the hopping term of the Ablowitz-Ladik model 



is 



j 



lim 

C — >oo 



*(0-C 



2Af 3 + l 



c 



2N.-1 



(42) 
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and 



i = lim 



t(0-C 



-2N S 



c 



1-2N. 



(43) 



whereas the interaction is related to the quantum deter- 
minant 



lndet,(T) = ^g(l+q}q j ) + (2N s + l)r 1 . (44) 

i=-N, 

In simple words: the hopping term (the kinetic energy) 
commutes with the interaction since the transfer ma- 
trix and the quantum determinant commute each other. 
Eigenvectors of the transfer matrix can be constructed 
by algebraic Bethe Ansatz (see Refs. 0, or Chapter 
VII of the book .3])- Let us denote the off-diagonal ele- 
ment of monodromy matrix (|35|) by C(£) = ?2i(C). The 
eigenvectors of monodromy matrix and Hamiltonian can 
be written as 



lx> = f[c(Cfc)|o) 



(45) 



fe=i 



Here |0 > is the Fock vacuum qj\0) = 0. The eigenvalue 
of monodromy matrix is 



_ /-2JV,+1 BT) 



o(0= c 



n 



sinh(A — A/j —77) 
k—1 sinh(A — Xk) 



1 /--2N s -l c nri TT™ sinh(A-A fc +?7) , A „s 
+ "» e 1 lfc=l s inh(A-A fc ) ' 

Here exp(A) = £ (the same for exp(Afc) = Cfe)- The vari- 
ables Afc have to satisfy the Bethe equations 



2(2N s +l)X a 



-A- sinh(A a — A fc + 77) 
JJ- sinh(A Q - Afc - ry) 



(47) 



The vectors (|45f) are also eigenvectors of the logarithm 
of the quantum determinant l|44|) with eigenvalues equal 
to (2Ns + l)r]. The energy levels [ eigenvalues of the 
Hamiltonianl are 



2(l-e 2 ")^cosh2Afc. 



(48) 



fc=i 



The operators qj can be represented in terms of the stan- 
dard Bose operators [a n , aj„] = 5^ 



'exp2r/(ata + 1) - 1 



y y 1 + ata 

The correlation functions of the model were evaluated in 
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